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We discuss the properties of the collective modes of a system composed by a thermalized 
quark-gluon plasma traversed by a relativistic jet of partons. The transport equations obeyed 
by the components of the plasma and of the jet are studied in the Vlasov approximation. 
Assuming that the partons in the jet can be described with a tsunami-like distribution 
function we derive the expressions of the dispersion law of the collective modes. Then the 
behavior of the unstable gauge modes of the system is analyzed for various values of the 
velocity of the jet, of the momentum of the collective modes and of the angle between these 
two quantities. We find that the most unstable modes are those with momentum orthogonal 
to the velocity of the jet, and the effect is stronger for ultrarelativistic jet velocities. Our 
results suggest a new possible collective mechanism for the description of the jet quenching 
phenomena in heavy ion collisions. 

§1. Introduction 

One of the methods for unveiling the properties of matter produced in ultrarela- 
tivistic heavy-ion collisions is to study the propagation properties of high px partons 
generated by hard scatterings in the initial stage of the collision. When the jet of 
partons travels across the medium it loses energy and degrades, mainly by radiative 
processes (see 1 - 1 for reviews). The energy and momentum of the jet are absorbed 
by the plasma and result in an increased production of soft hadrons in some given 
directions. 

We propose a novel mechanism 2 )' 3 ) for describing how the jet loses energy and 
momentum while travelling in a thermally equilibrated quark-gluon plasma (QGP). 
Since the jet of particles is not in thermal equilibrium with the QGP it perturbs 
and destabilizes the system inducing the generation of gauge fields, even if the jet of 
colored particles is neutral over a coarse grain volume. Some of these gauge modes 
are unstable and grow exponentially fast in time absorbing the kinetic energy of the 
jet. 

The novel mechanism we propose for the QGP is well-known in the context of 
traditional plasma physics. 4 ) The study of the interaction of a relativistic stream of 
particles with a plasma is a topic of interest in different fields of physics, ranging 
from inertial confinement fusion, astrophysics and cosmology. When the particles of 
the stream carry electric charge, plasma instabilities develop, leading to an initial 
stage of fast growth of the gauge fields. 

The study of chromo instabilities is a very active field of research nowadays (see 5 ) 
for a review), as it has been claimed that they may speed up the isotropization and 
equilibration process in an anisotropic system. Here we discuss another physical 
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setting where plasma instabilities may play also an important role in heavy ion 
collisions. 

In order to study the system composed by the jet and the plasma we have 
employed two different methods. In 2 ) both the plasma and the jet are described 
using a fluid approach. This approach developed in 6 ) has been derived from kinetic 
theory expanding the transport equations in moments of momenta and truncating 
the expansion at the second moment level. The system of equations is then closed 
with an equation of state relating pressure and energy density. The fluid approach 
has several advantages with respect to the underlying kinetic theory. The most 
remarkable one is that one has to deal with a set of equations much simpler than 
those of kinetic theory. Then one can easily generalize the fluid equations to deal with 
more complicated systems. This is a strategy that has been successfully followed in 
the study of different dynamical aspects of non-relativistic electromagnetic plasmas, 4 ) 
even if one is considering a system far away from its hydrodynamical regime. 

In the second approach we consider the same setting, i.e. an equilibrated plasma 
traversed by a relativistic jet, but we use kinetic theory instead of the fluid approach. 
Transport theory provides a well controlled framework for studying the properties of 
the quark-gluon plasma in the weak coupling regime, g <C 1. Indeed it is well known 
that the physics of long distance scales in an equilibrated weakly coupled QGP can 
be described within semiclassical transport equations. 7 )~ 9 ) In this approach the hard 
modes, with typical energy scales of order T, are treated as (quasi-)particles which 
propagate in the background of the soft modes, whose energies are equal or less 
than gT, which are treated as classical gauge fields. This program has been very 
successful for understanding some dynamical aspects of the soft gauge fields in an 
almost equilibrated QGP. 10 )' 11 ) 

In this talk we only review the results of the above second approach. Details 
of the comparison of the two approaches, that give similar results, can be found in 
Ref. 3 ) 

§2. Kinetic theory approach 

We consider a system composed by a quark-gluon plasma traversed by a jet 
of partons. We assume that the system is initially in a colorless and thermally 
equilibrated state and we will study the behavior of small deviation from equilibrium. 
The distribution function of quarks is 

Q(p,x) =f F q D (po) + SQ( P ,x), (2-1) 

where the various quantities are hermitian matrices in the fundamental representa- 
tion of SU(3) (we have suppressed color indices) and where f'p D (po) = ePi Jr +1 i s 
the (colorless) Fermi-Dirac equilibrium distribution functions. A similar decomposi- 
tion can be done for the distribution function of antiquarks Q(p,x), and for gluons 
G(p,x) that we will not write down here (see Ref. 3 ) for details). 

Also the jet is assumed to be initially colorless, but small color fluctuations are 
present: 

W iet (p, x) = f iet (p) + 6W iet (p, x) . (2-2) 
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For the initial jet distribution function we will not consider a thermal distribution 
function. We will indeed approximate the distribution function of the jet with a 
colorless tsunami-like form 12 ) 

fj et (p) = nu° 5^(p-Au) , (2-3) 

that describes a system of particles of constant density n, all moving with the same 
velocity = (u°,u) = 7(1, v), where 7 is the Lorentz factor and A fixes the scale 
of the energy of the particles. 

This distribution function is adequate for describing a uniform and sufficiently 
dilute system of particles, and it represents an extreme crude model of a jet of 
energetic partons, chosen to simplify our analytical estimates. Of course, it is possible 
to use more involved distribution functions, assuming that the density of particles 
composing the jet is not uniform and there is a spread in momentum. 

The distribution function of quarks satisfy the following transport equation 

pPDnQfax) + 9 -p» {F^(x),d»Q(p,x)} = C , (2-4) 

whereas the distribution functions of antiquarks, gluons and of the particles of the 
jet satisfy similar equations. Here g is the QCD coupling constant, with {...,...} 
we denote the anticommutator, dp is the four-momentum derivative, and is a 
covariant derivative. In Eq. (|2-4p C represents the collision term. However, for time 
scales shorter than the mean free path time the collision term can be neglected, as 
typically done in the so-called Vlasov approximation. 

The knowledge of the distribution function allows one to compute the associated 
color current, which in a self-consistent treatment enters as a source term in the Yang- 
Mills equation. It is important to note that at very short time scales the different 
components of the system formed by the plasma and the jet interact with each other 
only through the generated average gauge fields. 

The contribution to the polarization tensor of the particles of species a (where 
a refers to quarks, antiquarks, gluons or the partons of the jet) in the Vlasov ap- 
proximation is 

n ab,a( k ) = -9 C F 5 ab J faip) . ( 2 ' 5 ) 

where a, b are color indices and Cp is the value of the quadratic Casimir associated 
with the particle specie a which takes values 1/2 and 3 for the fundamental and 
adjoint representations of SU(S), respectively. The momenta measure is defined as 

/- = / 2®(Po)S(P 2 ~ ml ) , (2-6) 

where m a is the mass of the particle of specie a. For simplicity we assume that the 
particles belonging to the plasma are massless. Instead the particles of the jet have 
a non- vanishing mass. 



4 



Mannarelli and Manuel 



When f a (p) is a thermal equilibrated distribution function, Eq. (12-5P reduces to 
the form of the hard thermal loop (HTL) polarization tensor. 7 )~ 9 ** However for the 
tsunami-like distribution function, Eq. (|2-3p . the polarization tensor obviously takes 
a different form. 

The gauge fields obey the Yang-Mills equation 

D^(x) = 8f t (x) = 5f p (x) + SjV et (x) , (2-7) 
where we have defined 

= -§j/V [SQ( P ,X) - 6Q(p,x) + 2r a Tr[T a SG(p,x)]] , (2-8) 

which describes the plasma color current, and 

Sj^{x) = -~j^5W^{p,x), (2-9) 

which describes the fluctuations of the current associated with the jet. 

Equation (|2-7p together with Eq. (|2-4p form a set of equations that has to be 
solved self-consistently. Indeed the gauge fields which are solutions of the Yang- Mills 
equation enter into the transport equations of every particle species and, in turn, 
affect the evolution of the distribution functions. 

§3. Collective modes in the system composed by the QGP and jet 

We now consider the collective modes of the system composed by an equilibrated 
QGP traversed by a jet of particles. We are interested in very short time scales when 
the Vlasov approximation can be employed. The effect of the beam of particles is 
to induce a color current, which provides a contribution to the polarization tensor. 
The polarization tensor of the whole system is additive in this short time regime, 
meaning that 

nr(k) = n^(k) + n^(k), (3-i) 

where n p u (k) and nj*£(k) are the polarization tensor of the plasma and of the jet 
respectively. The total dielectric tensor is given by 



e?( Wj k) = + (3-2) 



and the dispersion laws of the collective modes of the whole system can be determined 
solving the equation 



det 



k 2 ^ -k l k? •-u 2 e\ j (k) =0. (3-3) 



The solutions of this equation depend on |k|, |v|, cos 6 = k • v, uj 2 = uj 2 + tv 2 et 
where uj p and a;j e t are the plasma frequency of the QGP and of the jet, respectively, 
and on 
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We find u? ei = g 2 fi/2A, so it depends on the density of the jet as well as the 
mean energy that their partons carry. 

Clearly, when the plasma and the jet do not interact, they have stable collective 
modes. However, once we consider the composed system of plasma and jet interacting 
via mean gauge field interactions, unstable gauge modes appear. 



0.03 




Fig. 1. Imaginary part of the dispersion law of the unstable longitudinal mode for the system 
composed by a plasma and a jet in the case k || v as a function of the momentum of the mode 
at b — 0.1 (left) and at 6 = 0.02 (right) for four different values of the velocity of the jet, |v|. 
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Fig. 2. Imaginary part of the dispersion law of the unstable mode for the system composed by a 
plasma and a jet with the kinetic theory approach for k _L v as a function of the momentum of 
the mode. Left panel refers to b = 0.1 and right panel refers to 6 = 0.02. In both case results 
for four different values of the velocity of the jet, |v|, are shown. 



In Figs. [T] and [2] we report the results for the unstable mode obtained with the 
kinetic theory approach for the case where k || v and k _L v, respectively, for different 
values of the jet velocity, and different values of b. 

From our results one concludes that the most unstable gauge field mode corre- 
sponds to the case k 1 v. The largest values of growth rates correspond to high 
values of jet velocities. One can also see that the growth rate increases when increas- 
ing the value of b. 

We can make an estimate for the time scale t for the development of the stream 
instabilities in the QGP. In the weak coupling limit, and for T ~ 350 MeV, we find 
t ~ 1 — 2 fm/c, meaning that it is a very fast phenomena. 
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§4. Energetic considerations 

With the generation of exponentially growing gauge fields it is clear that there 
must be an effective transfer of energy from the jet to the growing chromoelectro- 
magnetic fields. Such a transfer can be studied by analyzing the energy-momentum 
0^ v associated to both the partons of the jet and the plasma, as well as that of the 
gauge fields. 

The energy-momentum tensor due to the quasiparticles obeys 

Tr(^ t FH = ^r + ^fet (4-1) 

In collaboration with Michael Strickland 15 -* we are currently analyzing the dy- 
namical evolution of the energy of all the modes in the problem. For such a study we 
are employing the numerical code of Ref., 16 ) originally developed to study Weibel in- 
stabilities but now adapted to consider our physical setting. Our study is performed 
for a SU (2) non-Abelian group, and we use a Gaussian distribution to model the 
jet distribution function. The parametrization used for this function is such that 
in a given limit we reproduce the tsunami-like distribution function. The results 
of this numerical analysis will be presented elsewhere. 15 ) They confirm our initial 
expectations of jet energy loss. 

§5. Outlook 

Stream instabilities in the QGP can represent a collective mechanism for jet 
quenching. In a weak coupling scenario, this process is in principle much faster than 
collisional loss, which is suppressed by an additional factor of 1/g 2 . Still the strength 
of the effect we are discussing depends not only on g, but also on the ratio of densities 
of the jet and plasma, as well as in the jet velocity. It is not easy to predict how the 
mechanism we are proposing competes with collisional and radiative energy loss, so 
as to assess its relevance for heavy ion phenomenology. We leave this question for 
future studies. 

Let us stress again that our model of the jet of partons was extremely crude, but 
enough for preliminary estimates and also to study the dynamical evolution of these 
stream instabilities. One should localize the jet in space, with the probable associated 
effect of having a localized energy deposition. So far we have only considered the very 
short time evolution of the composed system. At sufficiently long times, collisions 
should be considered, by the addition of the corresponding terms in the transport 
equations. Collisions might stop the growth of the unstable gauge field modes, if 
those do not saturate due to non-linear non-Abelian interactions, and they represent 
a different sort of mechanism of jet energy loss. A complete analysis of jet quenching 
in the QGP should probably include all these different energy loss mechanisms. 
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